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ABSTRACT 

We report optical and X-ray observations of the high-field polar RX J1007. 5-2017 performed between 1990 and 2012. It has an orbital 
period of 208.60 min determined from the ellipsoidal modulation of the secondary star in an extended low state. The spectral flux of 
the dM3- secondary star yields a distance of 790 ± 105 pc. At low accretion levels, RX J1007. 5-2017 exhibits pronounced cyclotron 
emission lines. The second and third harmonic fall in the optical regime and yield a field strength in the accretion spot of 94 MG. 
The source is highly variable on a year-to-year basis and was encountered at visual magnitudes between V ~ 20 and V ~ 16. In the 
intermediate state of 1992 and 2000, the soft X-ray luminosity exceeds the sum of the luminosities of the cyclotron source, the hard 
X-ray source, and the accretion stream by an order of magnitude. An X-ray high state, corresponding to the brightest optical level, 
has apparently not been observed so far. 
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1. Introduction 

In polars, a subclass of cataclysmic variables (CVs), a mag- 
netic white dwarf accretes matter from its Roche-lobe filling 
low-mass main-sequence companion. The strong magnetic field 
channels the flow onto a small accretion spot on the surface of 
the white dwarf, where the energy not expanded in the accretion 
stream is released in form of hard X-rays, soft X-rays, and cy- 
clotron radiation. About 100 of these systems are known, most 
of them discovered as bright soft X-ray and EUV sources. The 
spectral characteristics of the spot emission depend sensitively 
on the magnetic field strength B in the ac cretion spot and on 
the density of the accreted matter (e.g., (e.g. lKuiipers & Pringld 
U982tlFischer & Beuermannl2001l) . The distribution of observed 
field strengths in polars ranges from 10MG to about 200 MG, 
with only three stars known to poss ess a field strength ex- 
ceeding 100 MG dSchwope et alj|2009l) . For comparison, single 
magnetic white dwarfs are known to possess field strengths up 
to 1000MG and display a higher incidence of very high field 
strengths. It is not known why such systems are lacking among 
the polars. Either post-common envelope binaries contain fewer 
white dwarfs with very high field strengths, or we fail to detect 
these systems once they reach the CV-stage. RX J1007. 5-2017 
(henceforth RXJ1007) is one of the few high-field polars. We 
present here a comprehensive analysis of its properties. 

RXJ1007 was discovered in 1990 as a bright, very soft X-ray 
source in the ROSAT All Sky Survey and was subsequently iden- 
tified wit h a 19-mag magnetic ca t aclysmic variable of the polar 
subclass dBeuermann & Thomas! 1 1 993b iThomas & Beuermann! 
ll998HReinsch et al.lll999TT In this paper, we summarize our ob- 
servations of RXJ1007 carried out over more than two decades, 
from its discovery in 1990 to 2012. The analysis of the early 
data suffered from the lack of a sufficiently accurate ephemeris, 
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preventing a proper phasing of the data taken over the years in 
different wavelength regimes. This problem has now been over- 
come and the entire set of so-far unpublished data is presented 
here. Our analysis is based on optical observations, comprising 
time-resolved spectrophotometry, spectropolarimetry, and pho- 
tometry. It includes an analysis of the ROSAT X-ray data and 
a re-anal ysis of the XMM-Newton X-ray data previously dis- 
cussed by iRamsav & Cropper] (120031) . 

2. Observations and data analysis 

Table Q] contains a log of our own and previously published ob- 
servations. RXJ1007 displayed substantial long-term variabil- 
ity both at X-ray and optical wavelengths. It was encountered 
in low states of accretion in 1997 and 2011/2012, in interme- 
diate states in 1992, 2001, and 2010, and in a high state in 
2000. The optical position of RXJ1007 is or 2 000 = 10 h 07 m 34^6, 
£2000 = -20°17'32". 

2. 1 . X-ray observations 

RXJ1007 was observed from 24-26 November 1990 in the 
ROSAT All-Sky-Survey (RASS) with the Position Sensitive 
Proportional Counter (PSPC) as detector for a total of 510 s 
(25 sightings). The highly variable source 1RXS J100734.4- 
201731 had a mean count rate of l.Octss -1 and a very soft 
spectrum. It was subsequently observed with the ROSAT PSPC 
on 17 November 1992 with a total exposure time of lOks and 
again foun d to be very soft with a lo wer mean count rate of 
0.37 cts s'. lRamsav & Cropped (T2003) observed the source on 7 
December 2001 with XMM-Newton with the European Photon 
Imaging Camera (EPIC pn) for 5.5 ks. This observation showed 
the source in a brighter state again and confirmed the very soft 
X-ray spectrum. 



1 



H.-C. Thomas et al.: The high-field polar RX J1007.6-2017 



Feb 2000 RXJ 1007-20 




4000 5000 6000 7000 

Wavelength [A] 



8000 



Fig.l. Mean spectra of RXJ 1007 in February 2000 (top), 
January 1992 (center), and March 1997 (bottom). 



2.2. Optical spectrophotometry and spectropolarimetry 

Phase-resolved spectrophotometry was performed in January 
1992, using the ESO/MPI 2.2-m telescope with EFOSC2 and 
either grism Gl (FWHM resolution 50 A) or grism G3 (FWHM 
resolution 10 A). Further low-resolution spectrophotometry was 
performed on 6 March 1997. On 3 February 2000, when the 
source was at the brightest level encountered so far, we per- 
formed circular spectropolarimetry, using the ESO 3.6-m tele- 
scope with EFOSC2 and a customer-supplied quarter-wave 
plate. All spectrophotometry was placed on an absolute flux 
scale using standard stars. Fig. Q] shows the mean spectra of 
RXJ 1007 in the 1992, 1997, and 2000 observations, which 
demonstrate its variability. In the 1997 low state, the red part of 
the spectrum is dominated by the TiO bands of the secondary 
star and the visual magnitude varies over the orbit between 
V = 19.4 - 20.0. The weak Balmer and the intense cyclotron 
emission lines indicate that low-level accretion is still taking 
place. The 1992 intermediate and 2000 high states are charac- 
terized by significantly increased levels of Balmer line emission 



Table 1. Summary of observations: X = X-ray, S = spectropho- 
tometry, C = circular spectropolarimetry, P = photometry, WL 
white light. 



Date 



Telescope/Instrument Range 



Resol. 



24-26 Nov 1990 
9-13 Jan 1992 

17 Nov 1992 
15 Nov 1995 

6 Mar 1997 
14 Dec 1999 

2-3 Feb 2000 

7 Dec 2001 
19-20 Jan 2010 
Jan -Mar 2011 

2 Jan 2012 



ROSAT/PSPC 
ESO 2.2m /EFOSC2 

ROSAT/PSPC 
ESO 2.2m/EFOSC2 
ESO 2.2m/EFOSC2 
2MASS 

ES0 3.6m/EFOSC2 
XMM/EPIC pn 
MONET/N 
MONET/N 
MONET/N 



X 0.1 1-2.0 keV 

S 3500-9000A 

S 3500-5400A 

X 0.1 1-2.0 keV 

S 3800-9 100A 

S 3500-1 0200A 

P J, H, K s 

C 3378-75 16A 

X 0.22-8.5 keV 

P IbcsscII, Sloan g 

P Wn,WL 

P iBessell 



40 A 

10 A 

30 A 
50 A 

15 A 



and the associated continuum, but only by a moderate increase 
in cyclotron emission, as discussed in Sect. 13.61 In the 2000 high 
state, the source varied over the orbit between V = 16 and 17. 

The profiles of the Balmer and helium emission lines in 
the medium-resolution spectra of 1992 and 2000 display phase- 
dependent asymmetries, but the narrow and broad line compo- 
nents, typical of polars, cannot reliably be separated. For sim- 
plicity, we measured radial velocities by fitting single Gaussians 
to the strong Balmer and Hell/14686 emission lines. Given the 
lack of better-resolved low-state spectra, we measured radial ve- 
locities and line fluxes also from the weak Ha emission lines 
as well as radial velocities from the unresolved near-infrared 
NaI/18190 absorption-line doublet in the 1997 spectra. In retro- 
spect, the phasing of the velocities suggests that the bulk of the 
intermediate and high-state emission originates from the mag- 
netically confined part of the accretion stream and the low-state 
Balmer emission from the heated face of the secondary star. 

2.3. Optical and IR photometry 

We performed time-resolved relative optical photometry in 
2010-2012, using the MONET/North telescope at McDonald 
Observatory via the MONET browser-based remote-observing 
interface. The data were taken with an Apogee ALTA E47+ 
lkxlk CCD camera mostly in Bessel / (central wavelength 
8000 A) and white light (clear filter) with exposure times 
of 60s. Photometry was performed relative to the compar- 
ison star USNOA 067510804546 (2000: RA=10 h 07 m 32.87 s , 
DEC=-20°1 5'36.3"), which has a DENIS /-band magnitude of 
13.14+0.03 dFouque et al.ll2000l) . Combined with standard-star 
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Eq. 1. 
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Fig. 2. Orbital variations of RXJ1007 from top: (1) 1990 X-ray 
light curve from the ROSAT All Sky Survey; (2) 1992 ROSAT 
pointed observation; (3) 1992 quasi [/-band flux; (4) 1997 
quasi t/-band flux; (5) 2000 visual flux; (6) 2000 Balmer/Hell 
emission-line radial velocity. The optical fluxes are in units of 
lO-^ergcnrV^ -1 . 



phase 

Fig. 3. Orbital variations continued: (1) 1997 Ho- emission-line 
flux; (2) 1997 NaLt8190 absorption-line radial velocity ; (3) 
1997 Ha emission-line radial velocity (open circles, dashed 
curve) and 1992 Balmer/Hell emission-line radial velocity 
(filled circles, solid curve); (4) 1992 / 750oA flux; (5) 1997/ 750oA 
flux; (6) 2010 /-band flux; (7) 201 1 /-band flux. 



measurements in the 1992 and 1997 campaigns, we derived a 
common /-band calibration for the spectrophotometry and pho- 
tometry. All our /-band data display a clear double-wave orbital 
modulation. RXJ1007 was found at / = 17.5 - 18.1 in January 



1992, /= 17.75-18.08 in March 1997, /= 17.28-17.74 in January 
2010, and / = 17.71 - 18.08 in January-March 2011 as well 
as January 2012. The source has an entry in the 2-micron All 
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Sky Survey (2MASS) catalog, which gives J= 16.30+0.10,//= 
15.63+0.12, and ^ s = 15.22±0.13. 



3. Results 

In this section, we describe first the effort to derive an accurate 
orbital ephemeris, which allows us to correctly phase the obser- 
vations collected over more than 20 years, and then continue to 
discuss the physical properties of the system. 

3. 1 . Orbital period 

Early in the project, we obtained approximate values of the or- 
bital period P of RXJ1007 by folding the X-ray fluxes (Fig. |2] 
panels 1 and 2 from top), the optical spectrophotometric fluxes, 
and the Balmer and Hell emission-line radial velocities (Fig. [2] 
panels 3-6, and Fig. [3] panel 3 from top, solid circles). The 
individual observations yielded P = 208.1 +0. 7 min (RASS), 
P = 205 ± 3 min (ROS AT pointed), and P = 208 .2 +0.3 min ( 1 992 
radial velocities), suggesting P^0.145d. All attempts to derive 
a long-term ephemeris from these data suffered from alias prob- 
lems. 

An alternative approach is suggested by the quasi-/ band 
fluxes in Fig.|3](panels 4 and 5 from top), which show what may 
be the ellipsoidal modulation of the secondary star in the 1992 
intermediate and the 1997 low-state data. The two light curves 
refer to the 7450-7550 A band, which measures the flux of the 
secondary star with some contribution from the white dwarf and 
the accretion stream. These early data were supplemented by the 
/-band light curves taken in the 2010 intermediate and the 201 1- 
2012 extended low states (Fig. [3] bottom two panels). For the 
dM3- secondary star in RXJ1007 (Sect. |33j» tne spectral flux at 
7500 A and the mean Bessell /-band flux agree within a few per- 
cent (Beuermann 2006), allowing an easy comparison. The four 
light curves shown have a common ordinate scale normalized to 
the peak flux in 2010, which corresponds to / = 17.3. Standard 
theory of the ellipsoidal modulation for Roche-lobe filling stars 
predicts two minima, of which the deeper one, according to von 
Zeipel's (1924) law, occurs at superior conjunction of the sec- 
ondary star. This order may be reversed if the low-gravity hemi- 
sphere of the secondary is radiatively heated by the accreting 
white dwarf. The different fluxes in the displayed light curves 
cannot be explained by heating alone, however, and the presence 
of an accretion-induced component is indicated for at least the 
1992 and 2010 intermediate states. 

In a first step toward an ephemeris, we corrected the ob- 
served minimum times in 1992, 1997, and 2010 for this addi- 
tional light source. We performed Fourier fits to the light curves, 
using two sinusoids with periods P/2 and P. The former approx- 
imates the ellipsoidal light curve (assuming two equal minima) 
and the latter represents the additional light source. The fit has 
five free parameters, the period P, the two amplitudes, and two 
phase shifts. Higher harmonics are not needed for the present 
purpose. The resulting fits are shown as solid lines in Fig. ^(bot- 
tom four panels) and an example of the P/2 sinusoid is added 
as a dotted curve in the third panel from the bottom. We used 
these fits for the sole purpose of correcting the observed min- 
imum times to those that the P/2 component alone would have 
had. In the low state, this correction practically vanishes. Table|2] 
lists the minimum times with this small correction included. In a 
second step, we performed a period search around P/2, using the 
times of Table [2] The best-fit period is f/2 = 0.07243 1961 (18) d 
with * 2 = 22.2 for 23 de grees of freedom (d.o.f.). The next-best 



alias corresponds to a difference of one cycle on the P/2 scale 
over one year. It has an unacceptable x 2 = 58.0 and can be ex- 
cluded. Hence, our fit fixes the orbital period P, but still leaves 
us with the choice of which of the two minima corresponds to 
inferior conjunction of the secondary. It is encouraging though 
that the deeper minima in 1992 and 2010 are separated by an 
even number of P/2 cycles and thus occur at the same orbital 
phase. We adopt the deeper minimum in the intermediate states 
as corresponding to inferior conjunction, a choice that we justify 
in Sect. 13.2.21 and obtain the ephemeris 

T = BJD(TDB) 2455215.96256(48)+0.144863 923 (36) E, (1) 

where Tq refers to the inferior conjunction of the secondary star. 
All orbital phases quoted in this paper refer to Eq. 1 . 

3.2. Orbital variations 

Figures [2] and [3] display the X-ray fluxes, optical fluxes, and ra- 
dial velocities folded over the orbital period. We determined the 
times of maximum flux or of maximum positive radial velocity 
by fitting sinusoids to the data, which are shown in the as dot- 
ted or solid curves in the appropriate panels. All time tags from 
the optical telescopes and from the ROSAT satellite are given 
in UTC, while XMM Newton has the leap seconds added and 
provides times in Terrestrial Time (TT). We converted all times 
to Barycentric Dynamical Time (TDB), which agrees with TT 
at the ms leveQ, and corrected them for the light travel time to 
the solar system barycenter, i.e., we use Barycentric Julian Days 
in TDB or BJD(TDB). The conversion from UTC was made 
with the interactive tool astroutil^J Note that we quote Julian 
days and not modified Julian days. 

3.2.1. X-ray flux 

The ROSAT X-ray light curves display an orbital maximum near 
phase <p — 0.78. A pronounced dip in the light curve of the 
pointed ROSAT observation (Fig. [2] second panel from top) oc- 
curs just before maximum, which lasts from </> = 0.73 to 0.78 and 
is likely produced by the accretion stream crossing the line of 
sight. Ramsay & Cropper (2003) observed RXJ1007 with the 
CCD cameras on board of XMM-Newton for part of the or- 
bital period at unknown phase. With the ephemeris of Eq.Q] we 
now find that their observation intervals covered orbital phases 
= 0.13-0.71 (EPIC MOS) and = 0.31 - 0.69 (EPIC pn), just 
missing the dip. The ROSAT light curves show X-ray emission at 
all orbital phases, suggesting either that the X-ray emitting spot 
never disappears completely behind the horizon, or that we are 
seeing emission from more than one region on the white dwarf. 
The X-ray flux varies substantially on time scales down to min- 
utes. Orbit-to-orbit variability is responsible for the low fluxes in 
the RASS light curve (Fig. [2] top panel, open circles). 

3.2.2. Radial velocities in the intermediate and high states 

Peak radial velocity of the Balmer and Hell/14686 line emission 
in the 1992 and 2000 intermediate and high states occur at orbital 
phases near X-ray maximum, suggesting that the line emission 
originates in the magnetically confined section of the accretion 
stream that leads from the stagnation point in the magnetosphere 
to the hot spot on the white dwarf (e.g. ISchwope et al.l 120051 

1 http://www.cv.nrao.edu/~rfisher/Ephemerides/times.html 

2 http://astroutils.astronomy.ohio-state.edu/time/ 
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IStaude et al.ll200ll) . The large mean radial-velocity amplitudes 
of the the H/3, Hy, H6, He, and He II A 4686 emission lines, 
A" =484 ±30 km s" 1 in 1992 and and K = 470 ±30 km s _1 in 2000, 
support this notion. At this point, it is appropriate to comment on 
the phase convention adopted in Section 13.11 Had we assigned 
<p — to the less deep /-band minimum, we would face the im- 
plausible situation that the radial velocity maxima would occur 
half an orbit offset from the X-ray maxima. Furthermore, the 
coincidence of maximum X-ray and 1992 and 1997 cyclotron 
fluxes would be destroyed. 

3.2.3. Optical-continuum fluxes 

The quasi-t/ band fluxes in the center two panels of Fig. ^rep- 
resent the cyclotron flux in the third harmonic, with a small 
contribution from the white dwarf in 1997 and an additional 
larger contribution from the accretion stream in 1992. Peak flux 
agrees in phase with the X-ray maximum and the radial veloc- 
ity maximum in 1992, and approximately also in 2000. We note 
that cyclotron beaming, which shapes the optical light curves in 
many polars, is not prominent for the low harmonics observed 
in RXJ1007. The quasi-V band flux / v i s observed in 2000 dis- 
plays an entirely different behavior, being in antiphase to the 
emission line radial velocities v ra d (bottom panels of Fig. [2j». The 
high-state continuum flux must, therefore, have an origin dif- 
ferent from the low- and intermediate-state cyclotron fluxes. A 
straightforward interpretation assigns the emission to the accre- 
tion stream as the dominant light source in the high state. The 
excess flux of the 2MASS data over the infrared continuum of 
the secondary star would then represent the Brackett and Pfund 
continua emitted in an intermediate or high state. The observed 
phase shift between f wls and v rac j in 2000 can be understood if 
the free-free and recombination continuum of the stream is op- 
tically thick and its brightest inner section is hidden from view 
at the phase of maximum recessional velocity. This model im- 
plies that cyclotron emission is not a major contributor to the 
optical flux in the high state. Cyclotron line emission, while the 
dominant source of the low-state U -band flux, is restricted to the 
optically thin fringes of the accretion spot in the high state and 
possibly shows little variation between the low and high states. 
This behavior is expected if high-density accreted matter carries 
its energy into subsphotospheric layers instead of expanding it in 
free - standing shocks above the photosphere of the white d warf 
(e.g. lKuiipers & Pringldll982l:lFischer & Beuermannll2001l) . 



Table 3. Times and phases of events around the orbit 



Maximum of 



BJD(TDB) 
24000000+ 



X-ray flux RASS 1990 48220.1588(40) 0.76(3) 

X-ray flux ROSAT 1992 48943.6190 (40) 0.83(3) 

Hfl-He, Hell radial velocity 1992 48631.7244(40) 0.81(3) 

Hfi-Hy, Hell radial velocity 2000 51577.6920(50) 0.91(4) 

Visual flux 2000 51577.6197(50) 0.41(4) 

Nal radial velocity 1997 505 13.7086 (70) 0.20(5) 

Htt radial velocity 1997 50513.7260(70) 0.32(5) 

Ha fine flux 1997 50513.7477(70) 0.47(5) 

Cyclotron flux 3670A 1992 48630.6984 (40) 0.73(3) 

Cyclotron flux 3750A 1997 50513.6496(30) 0.79(2) 
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Fig. 4. Observed spectrum in the 1997 low state at cyclotron 
minimum (green curve), fitted by a dM3- star and a blackbody 
(see text). 



3.2.4. Radial velocities and Ha line flux in the low state 

The top three panels in Fig.|3]display the Ho- line flux in arbitrary 
units, the NaI/18190 radial velocity, and the Ha radial velocity 
(open circles), as measured from the 1997 low -resolution spec- 
tra. The line flux peaks at = 0.5, when the irradiated face of the 
secondary is in view, and maximum recessional radial velocity 
occurs at <p 0.25, suggesting that these quantities trace the mo- 
tion of the secondary, providing additional support for our phase 
convention. Heating does occur, as indicated by the fact that the 
Nal line is best defined when the back side is in view. However, 
the Ha flux does not drop to zero at <p = 0, as expected for an 
origin on the heated side of the secondary and illustrated by the 
dashed curve calculated for an inclination i = 73° (Sect. 13.6b . 
Higher resolved spectra are needed to resolve the origin of the 
low-state Ha emission. The measured radial-velocity amplitudes 
are K^ Ha = 218 ± 30kms _1 and K' 2 NaI = 278 ± 30kms _1 . The 
true radial velocity amplitude Ko of the secondary is probably 
larger than the former value and smaller than the latter, and we 
compromise on A2 -250±20kms . 



3.3. The secondary star 

TiO features of the secondary star dominate the 1997 low-state 
spectra longward of 6200 A. The secondary star is best de- 
tected at orbital minimum of the cyclotron flux near = 0.25, 
which corresponds to the maximum of the ellipsoidal modu- 
lation. The mean of the three spectra at = 0.18, 0.26, and 
0.34 (Fig. |H green curve) shows the M star, the blue contin- 
uum of the white dwarf, and remnant cyclotron emission in the 
restricted intervals of 3700-4200 A (third harmonic) and 5500- 
6200 A (second harmonic). We fitted the observed spectrum in 
the intervals free of cyclotron emission, using representative 
spectra of M2 to M4 dwarfs from the Sloan Digital Sky Survey 
dRebassa-Mansergas et al.ll2.Q12b and a blackbody for the white 
dwarf. We find a spectral type dM3- with an uncertainty of half 
a spectral class (dotted black curve) and a blackbody tempera- 
ture of 22000 K (dashed black curve) for a white dwarf of 0.8 M Q 
with a radius of 7.4 x 10 8 cm at a distance of 790 pc (Sections l3.4l 
and !3.5l ). Extinction with A\ ~ 0.04, as suggested by the interstel- 
lar absorbing hydrogen column density derived from the X-ray 
spectral fits (Section [3^7] ). implies an /-band flux higher by 4%. 
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The 7500 A spectral flux of the secondary reduced to the level of 
the primary minimum at cf> — (Fig. [3] third panel from bottom) 
then is /7500 = 6.6 x 10 _17 ergcm _2 s _1 A _1 with an error of 10% 
based on the uncertainties in the flux calibration. Differential 
extinction would lead to a minimally earlier spectral type than 
adopted above. This value of /7500 is used in Sect. [33] to derive 
the distance ofRXJ1007. 

The amplitude of the ellipsoidal modulation is 25% of the 
peak flux, a number typical of a dM star viewed at high inclina- 
tion ;. For a dM3- star and i=73° (Sect. [3~6l l. the observed mod- 
ulation amplitude is reproduced using the appropriate gravity- 
darkening coeff icient B\ - 0.6 5 and limb-darkening coefficient 
u\ = 0.216 from lClarea d2000allbh . The ellipsoidal modulation is 
also seen in the wavelength interval A =4400-5200 A, where the 
secondary still contributes about half of the observed flux. 

3.4. Masses of the stellar components of RXJ1007 

The mean Roche-lobe radius of the secondary star is given by 
R2/R0 = 0.234 (M 2 /M ) 1/3 P 2/3 /, (2) 

where M 2 is its mass, P is the orbital period in hours, and 
/ - 0.989 (Kopal 1959). In CVs above the period gap, mass 
loss has usually driven the secondar y somewhat out of ther- 
mal equilibrium dKnigge et al.l 1201 lh . As a consequence, its 
radius is increased beyond that of an undisturbed field star 
by a higher perc entage than expected from tidal and centrifu- 
gal fo rces alone dKippenhahn & Thomas! fl970l: iRenvoize et al.l 
120021) . Depending on the long-term average of the mass transfer 
rate M, the secondary could have a mass as large as 0.40 M with 
a spect ral type dM2.9 or a s low as 0.24 M with a spectral type 
dM3.9 dKnigge et al.ll201ll their Table 2). The observed spectral 
type of dM3- favors a larger mass and settle on M 2 = 0.35 + 0.05, 
which allows for moderate bloating. From Eq. |2l we obtain 
R2 = (2.60 +0.13) 10 10 cm. An estimate of the primary mass Mi 
can be obtained from the radial-velocity amplitude K2 derived in 
Section[3X|] With / = 73°, we obtain Mi = 0.80+0. 15 M . 

3.5. The distance of RX J 1007 

We obtained the distance of RXJ1007 from the observed spec- 
tral flux of the secondary star combined with its radius and sur- 
face brightness. Roche geometry yields an effective radius of 
the projected cross section of the star at <p — and i = 73° 
of R pmi = 0.957 fi 2 = 2.49 x 10 10 cm. With the 7500 A spec- 
tral flux of the secondary star from above, /7500 = (6.6 + 0.7) x 
10~ 17 ergcni _2 s _1 A _1 , and the surface brightness of a dM3- star 
at750 0Aof F 7S qo = (6.4+D.9)xl0 5 ergcm- 2 s- 1 A- 1 (iBeuermannl 
I2006h . the distance is obtained as d — R VIO j (F 7500//7500) = 
(2.45+0.32) x 10 2I cm = 790±105pc. 

3.6. Cyclotron spectroscopy 

Cyclotron line emission dominates the 1997 low-state spectra 
and displays a marked orbital variation with a minimum at <p 
0.3 and a maximum at <f> =2 0.8. The latter coincides with the X- 
ray maximum (Fig. 0. Peak flux occurs when the line of sight 
is closest to the axis of the accretion funnel. This differs from 
the situation encountered in low-field polars, where cyclotron 
beaming results in a flux minimum in the direction of the funnel. 
In RXJ1007, however, we observed in low harmonics with little 
or no preference of emission perpendicular to the field (note that 
the first harmonic is emitted primarily along the field). 




Fig. 5. 2-D representation of the phase-dependent observed and 
calculated cyclotron lines of RXJ1007 in 1997. The best-fit cal- 
culated spectral flux (color) is compared with the observed fluxes 
(contour lines). The abscissa is frequency in units of 10 14 Hz, the 
ordinate is phase and data and model are displayed twice. The 
image shows the second and third harmonic with an indication 
of the first harmonic at the lowest frequencies. 



The phase-resolved cyclotron spectra were obtained by sub- 
tracting the contributions of the secondary and primary star from 
the 13 observed spectra of the 1997 data set. At peak flux, the 
second and third harmonics extend over the wavelength intervals 
5200-6700 A and 3500-4400 A, respectively. Between these in- 
tervals, the cyclotron flux is effectively zero. For a more detailed 
analysis, we rebinned the spectra to the spectral resolution of 
50 A and assigned 'errors' to the individual spectral bins that re- 
flect the scatter of the original data points within each 50 A bins. 
These spectra are displayed on a frequency scale in the 2-D im- 
age of Fig. |5] (contour lines). As an example, a single spectrum 
at maximum cyclotron flux is shown in Fig. [6] (top panel, green 
curve). 

We fitted the entire set of cyclotron spectra simultaneously, 
using model spectra calculated for an isot hermal homogeneous 
plasma slab dBarret & Chanmugamll 19851 so-called constant-A 
model). The free parameters of the multi-parameter fit were the 
magnetic field strength B, the plasma temperature T, the dimen- 
sionless slab thickness A, the solid angle O subtended by the 
accretion spot as seen from the Earth, the inclination i of the 
system, and the direction of the accreting field line given by the 
polar angle B and the the azimuth iff. The location of the spot 
is assumed to agree with the foot-point of the field line. We 
found that a consistent set of parameters exists that simultane- 
ously fits the thirteen spectra with 50 data points each. A color 



Table 4. Best-fit parameters for the fit shown in Figs.|5]and [6] 



Parameter 


Value 


Error 


Magnetic field strength (MG) 


94.0 


0.3 


Plasma temperature (keV)) 


2.1 


0.3 


Size Parameter A (10 4 ) 


2.8 


+2.3-1.3 


Solid angle n(l(T 26 sr) 


3.0 


0.3 


Inclination (°) 


73 


1 


Colatitude of field/spot B (°) 


12.3 


1.6 


Azimuth of field/spot i// (°) 


70 


4 
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Fig. 6. Representative cyclotron spectra showing the second and 
third harmonic in a field of 94 MG. Top: Cyclotron spectrum at 
<p = 0.8 in the 1997 low state (green) and theoretical cyclotron 
spectrum (black) obtained from a simultaneous fit to the entire 
set of spectra as shown in Fig. [3] Bottom: Circular polarization 
in the high state on 6/7 February 2000 measured around the flux 
minimum at orbital phase <p - 0. 



representation of the calculated spectral flux is shown in Fig. [5] 
The spectra represent the second and third harmonics in a field 
of 94 MG. The complete set of fitted parameters is provided in 
Table|4]along with their 1-cr errors. The fit has a x 2 = 5 17 for 643 
d.o.f. It is remarkable that a close-to-perfect fit to the observed 
data set is possible using the fairly simple theoretical model. The 
fit fixes the orbital inclination of RXJ1007 at i = 73° ± 1° (1-cr 
error), implying that the system just escapes eclipse of the white 
dwarf by the secondary star. With 2.1 keV, the plasma tempera- 
ture is lower than the mean temperature of 6 keV derived from 
the fit to the XMM spectrum, but both temperatures are substan- 
tially lower than expected for a free-standing shock on a white 
dwarf of 0.8 M @ . 

The dominance of emission from the magnetic pole with 
B = 94 MG persists from the low state into the high state of 
RXJ1007. Even against the intense Paschen continuum of the 
2000 high state, the second and third harmonics stand out in 
the circular polarization at practically the same wavelengths as 
in the low state (Fig. [6] bottom panel). The line shape, fur- 
thermore, suggests a low temperature similar to the 2keV of 
the low state, consistent with emission from the low-density 
fringes of the accretion spot, in which cyclotron emission is 
the prime cooling agent (Woelk & Beuermann 1996, Fischer 
& Beuermann 2001). Radiation-hydrodynamic calculations sug- 
gest that columns with different temperatures and plasma densi- 
ties may coexist in the spot and that cyclotron cooling dominates 
in low-density columns, in which the plasma cools before the ion 
and electron temperatures can equilibrate (Woelk & Beuermann 
1996, Fischer & Beuermann 2001). 

Despite the generally excellent agreement, some differences 
between the observed and calculated line profiles are notewor- 
thy. An observed feature that is not present in the blackbody- 
limited line flux is the pronounced central depression in the sec- 
ond harmonic around 6000 A. It is likely that this depression rep- 
resents the Fraunhofer-type absorption produced by the temper- 



ature stratification in colu mns that are optically thick in the emit- 
ted cyclotron harmonics dWoelk & Beuermannl l"l996). A similar 
profile was observed in the cyclo tron lines of the high-field po- 
lar UZ For dRousseau et alJl996h and was succ essfully modeled 
with the theory of Woelk & Beuermann! i 19961) . The slight notch 
in the long-wavelength fall-off of the model profiles results from 
differences in the line profiles for the ordinary and extraordinary 
rays in this fairly simple theory. 



3.7. The X-ray spectrum 

We have re-analyzed the mean X- ray spectrum taken with 
XMM-Newton on 7 December 2001 dRamsav & Cropperil20"oH 
their Fig. 8), using an improved spectral representation of the 
soft and hard X-ray compo nents. As demonstrated for AM 
Her dBeuermann et al]l2012l) . a model including two or more 
blackbody components may yield a more realistic estimate of 
the soft X-ray lu minosity than the single-bl ackbody assump- 
tion employed by iRamsav & Cropper! (120031) . Our fit involves 
two blackbodies with best-fit temperatures kTi = 59 + 3 eV 
and kTi = 32 ± 5 eV for a best-fit absorbing column density 
A^h = 1.9 x 10 20 H-atoms cm 2 . The hard X-ray component was 
modeled with two MEKAL spectra for solar abundances and 
temperatures kr = 6.0keV (fixed) and kr = 0.19 + 0.01)keV. 
The fit has a X 2 = U5 for 107 d.o.f. (x; = 1.07). The integrated 

leters are 3.5 x 10~" ergcm _2 s _1 
1 for the two blackbodies and 3.5 x 
1 for the two ther- 
mal components , respectively. Guided by the case of AM Her 
(IChristianl 120001) . we considered that the hard X-ray compo- 
nent of RXJ1007 possibly contains highly absorbed compo- 
nents that are not accounted for by the XMM spectral fit at 
E < 8keV. We estimated that this could contribute another 
(1 -2) x 10~ 13 ergcm _2 s _1 . Since the presence of a substantial 
hard X-ray component is uncertain, it is entered in Table [5] with 
a question mark. We caution, furthermore, that the correlated 
errors between the blackbody temperature and the interstellar 
absorbing column imply large uncertainties, preventing us from 
quoting errors for the fluxes. Despite all the uncertainties, how- 
ever, the emission of RXJ1007 is clearly dominated by soft X- 
rays. In the 1992 pointed ROSAT observation, the source was 
fainter by almost an order of magnitude. The fit parameters are 
less well-defined than those of the XMM fit and are not quoted 
here. 



fluxes for the best-fit parameters are 3.5 x 10 "ergcm "i 
and 8.5 x 10~ 13 ergcm~ 2 s 



10 13 ergcm 2 s 



1 and 8.5 x 10 13 erg cm V 



Table 5. Bolometric fluxes and luminosities of RXJ1007 for the 
intermediate states of 1992/2001. See text for the geometry fac- 
tors T]7l. 



Component 


F 




w 


L 




(erg cm~ 2 s 






(ergs" 1 ) 


Soft X-ray blackbody, 32 eV 


8.5 x 10" 


n 


In 


3.2 x 10 33 


59 eV 


3.5 x 10' 


1 1 


In 


1.3 x 10 33 


Thermal X-rays 0.2 keV 


9.0 x 10" 


13 


An 


6.7 x 10 31 


6.0 keV 


2.8 x 10" 


13 


\n 


2.1 x 10 31 


> lOkeV ? 


1.5 x 10' 


13 


\n 


1.1 x 10 31 


Cyclotron + accretion stream 


4.0 x 10" 


12 


An 


3.0 x 10 32 


Total 








4.9 x 10 33 
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Fig. 7. Spectral energy distribution of RXJ 1007-20. Left: IR-optical-UV part. Right: X-ray part. See text for further explanation. 



3.8. Overall spectral energy distribution 

Fig. [7] shows the spectral fluxes from the infrared to the hard 
X-ray regime collected into a single overall spectral energy dis- 
tribution (SED). For the chosen form v/ v vs. v, the integrated en- 
ergy flux of a component with spectral flux f v is F — J vf v dlog v, 
with log the natural logarithm. Given the flux of a spectral com- 
ponent, the luminosity L = rjnd 2 F was calculated with c/ = 790pc 
and the geomet ry factor 77 quoted in Tab le[5] As for the prototype 
polar AM Her (iBeuermann et al.ll2012h . we adopted emission of 
the soft X-ray blackbody and of cyclotron radiation into 2 ji and 
of the accretion stream and the hard X-rays into 4 n. 

The left panel shows our spectrophotometry of 1992, 1997, 
and 2000 (black curves), supplemented by the 2MASS J, H, and 
K-band fluxes of 1999 (red filled squares) and the visual and 
ultraviolet fluxes measured with the optical monitor on board 
of XMM-Newton simultaneous to the 2001 X-ray observation 
(red filled circles). For the low state, which is dominated by the 
seconda ry star, we added th e flux distribution of the dM3 star 
LHS58 dLeggettet al.lll996l green curve) adjusted to the 1997 
low-state spectrophotometry. The right panel shows the incident 
spectra for the 2001 XMM pn observation (red filled circles) 
and the 1992 ROSAT PSPC observation (open circles). For the 
XMM spectrum, we also included the 'source spectrum' cor- 
rected for interstellar absorption, using the X-ray spectral pa- 
rameters quoted above (dashed curve in the left and the right 
panel). 

There is a close agreement between our 1992 spectropho- 
tometry and the 2001 XMM optical-monitor fluxes in the vi- 
sual and ultraviolet bands, suggesting that they delineate the 
common SED of an intermediate state of accretion assumed by 
RXJ1007in 1992 and 2001. This SED is well documented from 
the infrared to the hard X-ray regime. About 92% of the bolo- 
metric flux in this state is emitted as soft X-rays, 6% as cy- 
clotron radiation and stream emission, largely in the ultravio- 
let, and only 2% as hard X-rays. The dominance of soft X-rays 
implies that most of the accretion energy is released in shocks, 
which are buried in the photosphere of the white in a scenario 
first described bv lKuiipers & Pringldd 19821) . The luminosity ra- 



tio L softX /(L cyc + Lhardx + Stream) ~ 11. The low-level cyclotron 
emission in 1997 suggests that the associated X-ray flux was low 
as well. Similarly, we can only speculate that the optical high 
state in February 2000 was accompanied by a correspondingly 
increased soft X-ray flux. The rough equality of the ROSAT and 
XMM X-ray fluxes in the l-2keV range hints at a lower vari- 
ability of the hard X-ray fluxes, consistent with the notion of 
iKuiipers & Pringlel dl982h that an increased M arises primarily 
from dense blobs, which penetrate to subphotospheric layers. 

For a white dwarf with Mi = 0.8 M and Ri = 7.4 x 10 s cm, 
the accretion luminosity in the 1992/2001 intermediate state of 
L^c — 4.9 x 10 33 ergem^s -1 requires an accretion rate M = 
RiU c J(GMi) 3.4 x lO^gs- 1 or 5.4 x 10 -10 M yr _1 . The 
lack of X-ray coverage prevents an estimate of the accretion rate 
in a high state, but if it follows the optical continuum, it could 
exceed 10~ 9 M y r 1 and, thereby, rea ch into the realm typical of 
long-period CVs dKnigge et al.ll201 ll) . 

4. Conclusions 

We have presented a comprehensive study of the high-field po- 
lar RXJ 1007, which includes observations performed over more 
than 20 years. It is a highly variable source that was encountered 
at optical flux levels differing by 4 mag over the years. In the 
intermediate state of 1992 and 2001, RXJ 1007 has a soft X-ray 
luminosity L so ftx ~ 0.90 L acc 1 .2 L Q . In a high state of accretion, 
the soft X-ray luminosity is possibly higher still, both relatively 
and absolutely. With a polar field strength of 94 MG, RXJ 1007 
belongs to the few high-field polars. The main accreting pole 
is permanently in view and is the source of the cyclotron lines 
observed at all accretion levels. The observed cyclotron line pro- 
files and the low temperature in the emission region are perfectly 
consistent with cyclotron theory. 
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RXJ 1007 collected in early January 2012, before he suddenly died of heart fail- 
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with stellar structure theory and with close binaries. Being one of the earliest 
collaborators of Rudolf Kippenhahn, Hans-Christoph Thomas was the first to 



8 



H.-C. Thomas et al.: The high-field polar RX J1007.6-2017 



follow the evolution of a solar-like star numerically through the helium flash. 
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